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ABSTRACT 

We use ~12,000 spectra of ~3,500 magnetically active M0-M9 dwarfs from the Sloan Digital Sky 
Survey taken at 10-15 minute intervals, together with ~300 spectra of ^^60 M0-M8 stars obtained 
hourly with the Hydra multi-object spectrometer, to probe Ha variability on timescales of minutes to 
weeks. With multiple observations for every star examined, we are able to characterize fluctuations in 
Ha emission as a function of activity strength and spectral type. Stars with greater magnetic activity 
(as quantified by Lna/Lboi) are found to be less variable at all spectral types. We attribute this result 
to the stronger level of persistent emission in the high activity stars, requiring a larger heating event 
in order to produce measurable variability. We also construct Ha structure functions to constrain 
the timescale of variability. The more active objects with lower variability exhibit a characteristic 
timescale longer than an hour, likely due to larger, longer lasting heating events, while the less active 
objects with higher variability have a characteristic timescale shorter than 15 minutes. 

Subject headings: stars: activity — stars: chromospheres — stars: flare — stars: low-mass — stars: 
magnetic fields 



1. INTRODUCTION 

M dwarfs, the least massive and most numerous 
stars in our Galaxy, provide remote laboratories for 
observing high-energy magnetic phenomena in stellar 
atmospheres. The red and dim nature of M dwarfs 
makes the Ha spectral feature the most accessible tracer 
of magnetic activity in their convective outer la yers 
(|Bonn fc Schmitd [19781 : IWalkowicz fc HawlevI 120091) . M 
dwarfs with Ha in emission are classified as magnet- 
ically active, with a greater fraction of stars showing 
Ha emission at l a ter spe ctral types (Hawlcv ct al. 1993 
IGizis et all 120001: iWest et al.ll201lO . Magnetic activity 
is correlated w ith energetic flaring events tha t last sec- 
onds to hours (|Moffetti 11971 iKowalski et"all 120101 ) and 



photometric brightnes s increases as la.rge a s 8.2 magni- 
tudes in the i3-band (jAlmeida et al.l 120111) . The inci- 
dence of Ha emission and large flares decreases for stars 
farther from th e Galact ic Plane (West et al..,2006. 2008; 
Kowalski et al. 2009; Hi lton et al. l l2010l:IWalkowiCT et al.l 



20111) . This is interpreted as an age effect, leading to 
activity timescales of 1-8 billion years for early-late M 
dwarfs respectively. 

Variability in quiescent (non-flaring) Ha emission has 
been attribut ed to both stellar rotation and M dwarf ac- 
tivity cycles. iBopp k, Schmitzl (|l978l) related sequential 
Ha measurements to orbital phase and suggested that 
the observed changes were caused by the rotation of lo- 
calized active regions in and out of view on timescales of 
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days. Activity cycles with per iods of a few years have 
been ide ntified for M dwarfs bv Cinc unegui et al.l (|2007D 
and iBuccino et al.l (2011), and may affect the observed 
strength of Ha emission . Although l arge fiares are rela- 
tively infrequent events HHiltonll201lD , sporadic low-level 
magnetic field reconnection events ( "microflares" ) may 
also contribute to inherent Ha emission variations in ac- 
tive M dwarfs. 

There have been some recent investigations of Ha vari- 
ability on timescales of minutes to days. IGizis et "all 
(|2002l) measured the Ha equivalent widths (EWHA) 
for a sample of more than 600 active M0-M6 dwarfs 
from the Palomar /Michi gan State Univer sity (PMSU) 
survey of nearby stars ()Reid et "all '1995), and found 
that the variation in repeated measurements was near 
15% for stars with mean EWHA, (EWHA) < 5A, while 
the variability for stars with (EWHA) > 5A was in 
excess of 30%. In their study of EWHA variations 
on later type stars, iLee et ahl (|2010( ) identified signif- 
icant variability in ~80% of their sample of 43 M4- 
M8 stars. They found that the majority of variabil- 
ity events had timescales of 30 minutes or longer, and 
that the later-type stars were more variable, as mea- 
sured by i?(EW) = max(EWHA)/min(EWHA). Simi- 
larly, iKruse et al.l (|2010[ ) found an increase in the same 
metric (which they refer to as Rew) at later spectral 
types for timescales from 15 mi nutes to an hour in the 
Sloan Digital Sk y Survev (SDSS: lYork et al.l l2Ci00l) Data 
Release 7 fDR7: lAbazaiian et al.ll2009 t) component spec- 
tra. They also noted that Rew was larger for intermit- 
tent sources. 

Our analysis uses a set of spectra with exposure times 
of '^1 hour from the multi-object Hydra spectrograph on 
the WIYN telescope at Kitt Peak National Observatory, 
as well as the component spectra from SDSS, to extend 
the evaluation of Ha variability using new metrics and 
a broader range of timescales. We employ a construct 
known as the structure function, which is commonly used 
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Figure 1. Sample high-S/N active M4 dwarf spectra from the 
Hydra (top) and SDSS (bottom) data are shown. The data have 
moderate spectral resolution (2000-4000) and wavelength coverage 
from ~4000-9000A for SDSS and --4000-7000A for Hydra. 



in quasar variability analyses (e.g. iSimonetti etaTlllQSSl : 
iMacLeod et al.|[2011[ ). to further constrain the timescales 
of Ha emission yariability. 

The spectroscopic M dwarf samples we used are de- 
tailed in Section 2, along with our methods for determin- 
ing spectral types and measuring Ha emission strength. 
Analysis of Ha emission variability with spectral type, 
activity strength and timescale is presented in Section 3, 
and we summarize our findings in Section 4. 

2. DATA 

Our data come from two sources that are discussed in 
detail below: the Hydra instrument on the WIYN tele- 
scope (Section 2.1), and the Sloan Digital Sky Survey 
(Section 2.2). Representative spectra from each sample 
are displayed in Figure [TJ The spectral type distributions 
of the active M dwarfs in the two samples are shown in 
Figure [5] These moderate-resolution data sources (Hy- 
dra: R-4000, SDSS: R-2000) are useful for our Ha vari- 
ability analysis because they each provide multiple spec- 
tra of the same objects at well-defined time intervals. 

2.1. WIYN Hydra Data 

Our analysis employs data that were taken over 27 
nights between January 2002 and July 2008 as part 
of the Chandra Multiwavele ngth Plane Survey (ChaM- 
Plane: iGrindlav et al.l [20051 ) . which was designed to in- 
vestigate the fraction of diffuse X-ray emission in the 
Galactic plane due to unresolved X-ray point sources. 
Candidate point sources were identified using VRIHa 
photometry obtained with the Mosaic cameras at the 
Cerro Tololo Inter- American Observatory and the Kitt 
Peak National Observatory 4-meter telescopes. Spectra 
were obtained with the Hydra multi-object fiber-fed spec- 
trograph on the 3.5-meter WIYN telescope for objects 
with Ha — R < —0.3, indicating excess Ha emission, 
as well as optically detected objects with corresponding 
X-ray emission in archival Chandra images. Although 
the survey was designed to identify accretion powered 
sources such as cataclysmic variables and low-mass X- 
ray binaries, initi al analysis by the ChaMPlane team 
(jRogel et al.l I200d) revealed a sizable population of M 
dwarfs in the sample, whose spectra and photometry we 
utilize in this study. 

We assigned spectral t ypes to the Hydr a data us- 
ing the Hammer software ()Covev et al.l 120071) in manual 
mode, because the flux calibration and signal-to-noise 
ratio (typically ■^lO) of the spectra, especially at blue 
wavelengths, did not allow automatic typing. Each ob- 
ject was vis ually compared to the H ammer set of M dwarf 
templates ([Bochanski et al.ll2007D . The resulting spec- 
tral types are accurate to within 1-2 types. We verified 
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Figure 2. The distribution of spectral types is similar for the 
magnetically active M dwarfs in the SDSS (solid line) and Hydra 
(dashed line) samples. 

that the spectral types were broadly consistent with the 
measured V — I colors of the objects. 

The strength of Ha emission in each spectrum was 
quantified by measuring the equivalent wi dth of the 
line. W e followed the procedure described in lWest et al.l 
(2004') and adopted their definitions for the Ha line and 
continuum regions. Figure [3] identifies these regions for 
the four Hydra spectra obtained for a representative M4 
dwarf. We visually inspected every spectrum in the sam- 
ple to ensure that the EWHA measurements were not 
contaminated by cosmic rays or other artificial defects. 

Because the equivalent width is measured with respect 
to the local continuum level, it can be misleading to com- 
pare EWHA among stars of different spectral type, and 
hence continuum flux near Ha. We therefore computed 
the ratio of the Ha line luminosity to the bolometric lu- 
minosity (LHr t/Lhoi) for each object using the x conver- 
sion factors of lWest fc Ha wlev (2008*). This quantity rep- 
resents the continuum-independent activity level, which 
allows for Ha emission to be compared between stars of 
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Figure 3. Ha line profiles in the Hydra data for repeated ob- 
servations of an M4 dwarf. The solid vertical lines indicate the 
boundaries of the line region used in the equivalent width calcu- 
lations, and the dashed lines mark the pseudo-continuum regions 
on either side. The measured equivalent width for each spectrum 
is indicated to the left of the emission line, and the time elapsed 
since the beginning of the first exposure is shown on the right. 



Ha Emission Variability in Active M Dwarfs 



3 



different spectral types. 

Our analysis was limited to the chromosp herically ac- 
tive M dwarfs, identified by the criteria of iWest et"al] 
(|2011| ) that define detectable Ha emission. Because we 
intend to use these data in a time-domain analysis, we 
further restricted the sample to include only those stars 
with at least three spectra that were identified as active. 
Our final Hydra sample contains 312 spectra from 62 ac- 
tive M0-M8 dwarfs, and allows us to probe Ha emission 
variability on timescales from one to three hours. 

2.2. SDSS Spectra 

All SDSS spectra are composites of several (typically 
3-5) spectra with ^10-15 minute exposure times. The 
individual component spectra (with typical signal-to- 
noise ratio ^7) were initially made public as part of the 
SPSS Data Release 6 (DR6; lAdelman-McCarthv etHI 
120081 ) . Previous studies have examined the spectral vari- 
ability of M dwarfs usiiig the individual S PSS spectra 
from PR6 ()Hilton et all [20Tol) and PR7 (jKruse et al.l 
our analysis ex pands on these efforts. We use 
the IMlton et al.l (|2010( ) sample, which has the spectral 
type for each ob iect assigne d from the composite spec- 
tra presented in I West et all (j2008| ). The EWHA mea- 
surements for the component spectra followed the same 
prescription as for the Hydra data. Lna/Lboi values were 
obtained from the EWHA calculated for each object, and 
we verified that they agreed with those measured from 
the composite spectra given in IWest et al.l (|2008!) . We 
further restricted our sample to include only the persis- 
tently active objects that showed measurable Ha emis- 
sion in ev ery individua l exposure. We note that this re- 
duced the lHilton et al.l (2010) active sample by less than 
3%. 

Our final SPSS sample consists of 11,989 spectra for 
3,469 active M0-M9 dwarfs, and allows us to examine 
variability on timescales from ~15-100 minutes. Nearly 
700 stars were observed on multiple days, allowing us 
to further probe Ha emission changes on timescales of 
~l-20 days. 

3. VARIABILITY ANALYSIS AND RESULTS 

We first consider Ha variability as a function of spec- 
tral type, comparing our results with variability indica- 
tors used in previous analyses (Section 3.1). Our pre- 
ferred metric for variability is ctewha/ (EWHA), and we 
investigate this metric as a function of activity strength 
(Lna/Lboi; Section 3.2). Finally, we compute structure 
functions to investigate characteristic timescales in the 
low- and high-variability subsets of our data (Section 
3.3). 

3.1. Ha Variability with Spectral Type 

iLee et al.l (pOlOl ) and iKruse et al.l (pOlOl ) both found 
that Ha variability increases at later spectral type, using 
the Rew metric. Figure S] illustrates Rew as a func- 
tion of spectral type for our Hydra and SPSS samples. 
In agreement with their analyses, we find that the dis- 
tributions are not symmetric, but have significant exten- 
sions to large values of Rew, with these extensions being 
more pronounced at later spectral type. However, in our 
samples the median values (the median is preferred as 
a statistic instead of the mean for these skewed distri- 
butions) scatter around a roughly constant value. This 
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Figure 4. The Rew variability metric is shown as a function of 
spectral type for the SDSS data (light gray dots) and Hydra data 
(light blue diamonds), illustrating distributions with large numbers 
of low- variability objects and extended tails toward higher variabil- 
ity that are more noticeable at later types. The distributions are 
artificially broadened at each spectral type to more clearly show the 
actual number of objects. The median of each type is given by the 
black asterisks (SDSS) and blue filled circles (Hydra). There is no 
obvious trend in the medians with spectr al type. The r ed sq uares 
and purple triangles denoting results from IKruse et aLl 1 120101 ) and 
[Lee et al. (2010) respectively, are shown for comparison. See text 
for discussion. 



is in c ontrast to the IKruse et al.l ()2010[ ) and ILee et ahl 
(|2010D results, shown as red squares and purple triangles 
respectively. We attribute the difference to t he cri teria 
used to id e ntify the active stars. IKruse et al.l (|2010( ) and 
iLee et aP ()2010f ) include stars that they characterize as 
"intermittently" variable (i.e. only showing Ha emission 
in a subset of their spectra), and they exclude weakly 
active stars, which they classify as inactive. These dif- 
fer from our analysis as a result of their non-uniform 
equivalent width measurement criteria. Their intermit- 
tently variable stars are more prevalent at later types and 
systematically exhibit larger values of the Rew metric. 
We are able to confidently identify Ha at much smaller 
equivalent width, even in late-type M dwarfs, provided 
the spectra have sufficiently high signal-to- noise ratios 
(jWest et al.ll2008[) . Evidently, including these less ac- 
tive, but persistent, objects leads to the approximately 
flat distribution outlined by the SPSS and Hydra medi- 
ans in Figure 21 

3.2. Ha Variability with Activity Strength 

'Giz is et "ahl ([2002') used the standard deviation of re- 
peated EWHA measurements {aswHA) for each star as 
a metric for line strength variability, and found that this 
value increased with (EWHA) for their sample. We in- 
vestigate that result in Figure [S^, which illustrates the 
variation of aEWHA with (EWHA) for stars in both the 
Hydra and SPSS samples. The median aswHA for bins 
containing roughly equal numbers of stars are also shown. 
The abrupt cutoff at (EWHA) ~ lA is due to our re- 
quire ments for measur able Ha equivalent width in these 
data (|West et al.ll2008( ). 

We see that the median aEWHA is near zero for low 
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Figure 5. a) The standard deviation in EWHA (uewha) ver- 
sus mean EWHA ((EWHA)) is shown for each star in the Hydra 
and SDSS samples. Medians of representative bins are also shown 
(black asterisks for SDSS, blue circles for Hydra). There is a rising 
trend toward larger a and higher scatter as (EWHA) increases, b) 
The medians of the binned data for the fractional variability met- 
ric (Tgj4^£fx/(EWHA) are shown grouped by spectral type for the 
SDSS data (M0-M2 as green squares, M3-M5 as purple diamonds 
and M6-M9 as red crosses). There is a clear separation by spectral 
type, and a general decreasing trend toward lower variability at 
larger (EWHA). The median SDSS and Hydra data have the same 
symbols as in the top panel. 

(EWHA), rising to -O.SA at (EWHA) ~ 3-6A and 
increasing to '-I.SA at (EWHA) > 10 A. Our well- 
sampled data show a smooth, continuously increasing 
scatter of (Jewha with (EWHA) for the median points, 
in qualitative agreement with the much smaller sample 
of iGizis et all (|2000D . We also note that the maximum 
value of a EWHA increases with (EWHA). 

As noted above, the absolute equivalent width values 
depend on spectral type due to the changing continuum 
flux near Ha, with earlier type stars having relatively 
smaller equivalent width due to their higher continuum 
flux. We illustrate this effect in Figure [5}d using the frac- 
tional variability metric ctbivha/ (EWHA). The skewed 
nature of the a distributions at a given (EWHA) , evident 
in Figure [5^ (and analogous to the skewed distributions 
in Rew shown in Figure S]), again leads us to use the 
median value to characterize a ew h A / (ESSTRK) for each 
(EWHA) bin. 

The spectral type dependence of (EWHA) is clearly 
indicated by the separation of the early (M0-M2, green 
squares), mid (M3~M5, purple diamonds) and late (M6- 
M9, red crosses) types for the SDSS sample. The me- 
dian a Ew H A I (ES^^Pi) values for all spectral types are 
shown as black asterisks (SDSS) and blue filled circles 
(Hydra). The variability metric decreases with increas- 
ing (EWHA) from -0.2 for stars with low EWHA to 
~0.14 for stars with (EWHA) > SA. 
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Figure 6. The same medians of the fractional variability metric 
(o'£;M''Ha/(EWHA)) from Figure [sJj are shown as a function of 
activity strength log ^Hal^bol- The symbols are the same as in 
Figure [5}3. The downward trend in the SDSS data shows that 
variability decreases in more active M dwarfs, regardless of spectral 
type. The Hydra data do not show a coherent trend, which we 
attribute to poor sampling and a bias toward very active, early 
type stars. 

A better measure of magnetic activity strength that 
accounts for the continuum dependence on spectral type 
is Lua/Lboi- Figure [5] illustrates the same fractional vari- 
ability metric a ew h A / (J^^^^) as in Figure [5}), versus 
activity strength. The spectral types are no longer sepa- 
rated; in fact the three spectral type groups from Figure 
[5}d, have shifted horizontally so that they now overlap 
and form a smooth distribution. It is clear that the un- 
derlying variability relationship is not a function of spec- 
tral type but rather of activity strength. However, the 
median fractional variability (black asterisks) decreases 
with increasing activity from —0.3 at log(LHQ/L;,o;) — 
—4.5 to —0.1 at log(LHa/Lboi) — —3.5. The median value 
for the variability metric across the entire sample is 0.16, 
shown as the dashed line on the Figure. 

The Hydra data (blue filled circles) in Figures [5)d and 
[6] show relatively flat distributions that are influenced 
by small numbers, uneven spectral type sampling, and a 
bias toward very active stars especially at early spectral 
types (due to the Ha — R color selection of the targets) . 
Therefore, they are not as useful as the SDSS data for 
understanding the behavior of the variability metric with 
activity strength. 

The picture we have developed is that highly active 
stars of all spectral types are relatively less variable than 
low activity stars. This only becomes evident when 
the well-known underlying relationship between activity 
strength and spectral type is removed, as in Figure IHl 
Physically, the strong level of persistent emission in the 
high activity stars may be due to a large surface cover- 
age of active regions, while the low activity stars may 
have a relatively smaller filling factor. An increase in 
Ha emission, for example due to a new active region ap- 
pearing on the surface, or a magnetic heating episode or 
microflare in one of the existing active regions, will then 
manifest as barely visible in the more active star, while it 
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Figure 7. Three time bins from the structure function of the low-variabihty SDSS sources. Each point on the structure function in Figure 
[8] is a normahzed RMS (denoted by the vertical dashed lines) calculated from a distribution of AEWHA/(EWHA) that is populated by 
pairs of observations separated by a specific range of time, At (given in hours at the top of each panel). 



will contribute significant new emission in the less active 
star. This would naturally explain why the less active 
stars are relatively more variable. 

3.3. Structure Functions and Timescales of Ha 
Variability 

Structure functions are often used to explore the vari- 
ability timescales present in time-resol v ed observations 
of qu asars (e.g. iSimonetti et al.l 119851 : iMacLeod et al.l 
1201 It) . They represent the change observed between two 
measurements of some variable quantity, typically broad- 
band flux in the case of quasars, as a function of the 
time between measurements (At). It is expected that 
if the observed variability operates on a characteristic 
timescale (r), then for At < r, the measurements will be 
correlated. The variability amplitude will increase with 
greater time separation as the observations sample more 
of the total variability range until At = r. For At > r, 
the structure function maintains a constant amplitude as 
the measurements are uncorrelated and randomly sample 
the full variability range. 

We employ structure functions here to quantify the 
typical change in the equivalent width measurements for 
active M dwarf spectra in many bins of time separation, 
At, and to constrain the characteristic timescale of vari- 
ability, r. We divide our samples at the median value 
of the fractional variability metric in Figure |6] into the 
less active but more- variable stars {aEWHA/ (EWH-^) 
> 0.16), and the more active but less-variable stars 
{aEWHA/CEWRA) < 0.16). We then calculate the 
change in Ha emission as a fraction of the mean value, 
AEWHA/ (EWHA), between every pair of repeated mea- 
surements for the same M dwarf. 

To illustrate how the structure function amplitude is 
calculated, the distributions for three SDSS time bins 
from the less-variable (more active) subsample are pre- 
sented in Figure [T] Each time bin (the At range in- 
dicated on the Figure) contributes one point on the 
structure function. The typical change between mea- 
surements for a given At bin is taken as the RMS of 
the AEWHA/ (EWHA) distribution, normalized by the 
number of samples in that time bin. This accounts for 
the uneven sampling of time separations in our datasets. 



The normalized RMS values for these three bins increase 
from 0.13 to 0.16 with increasing At. The SDSS struc- 
ture functions were constructed from similar normalized 
RMS values that were calculated for 14 time bins span- 
ning the entire range of At sampled by our data. The 
time bins were chosen to provide good statistics at short 
At intervals (15 minutes to 2 hours), where there are 
many observations, and to cluster around typical At val- 
ues for longer intervals {e.g. one day, several days). The 
Hydra structure functions were only computed for a few 
time intervals, reflecting the approximately hourly ca- 
dence of the observations. 

Figure [5^ illustrates the SDSS (black asterisks) and 
Hydra (blue filled circles) structure functions for the 
high-variability sample {ctewha/ (E^H^) > 0.16), and 
Figure [8b gives the structure functions for the low- 
variability sample {a ew h A / (EWRA) < 0.16). The 
timescales range from 15 minutes to 20 days. The irregu- 
lar binning in At reflects the non-uniform separations of 
the measurements in time due to the sampling cadences 
of the surveys. We also require that the shortest time 
bins (At < 2 hours) in the SDSS structure function con- 
tain at least 150 measurement pairs, in order to obtain 
well-defined normalized RMS values (see Figure ^ . At 
larger time separations, the sample sizes are smaller and 
the error bars on the structure function amplitudes are 
correspondingly larger. 

The low-variability structure function for the SDSS 
sample in Figure |8b shows a significant trend with time, 
indicating increasing variability up to a timescale of at 
least an hour. The plateau at timescales shorter than 
At ~ 0.4 hours suggests that measurement uncertain- 
ties may dominate any variability amplitude effects be- 
low that time separation. 

The high- variability structure function (Figure UK) , on 
the other hand, shows a flat distribution at all times for 
the SDSS sample. The structure function amplitude is 
much larger, which is indicative of the higher level of 
variability seen in this sample. The measurement un- 
certainties are much less than the detected RMS (struc- 
ture function amplitude), indicating that the variability 
timescale is shorter than our smallest time separation 
bin, ~15 minutes. 
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Figure 8. The structure functions for a) the high-variability (less 
active) stars; and b) the low-variability (more active) stars, for the 
SDSS (black asterisks) and Hydra (blue filled circles) samples. The 
low-variability stars show an increasing structure function ampli- 
tude between 0.4-1 hour, indicating a timescale of more than an 
hour for the variability. The high-variability stars apparently vary 
on timescales shorter than our smallest time separation (~15 min- 
utes). The Hydra data are discussed separately in the text. 

While sparsely sampled, the Hydra data reinforce this 
interpretation. In Figure [8^, both the SDSS and Hydra 
exposure times are longer than the underlying character- 
istic timescale for the high-variability data, and there- 
fore all measurement pairs are randomly sampling the 
intrinsic equivalent width variations, leading to flat struc- 
ture function amplitudes. The longer exposure time for 
the Hydra data (an hour compared to 10-15 minutes for 
SDSS) damps the stochastic variations between expo- 
sures, and the amplitude of the Hydra structure func- 
tion is correspondingly smaller. The low- variability Hy- 
dra sample in Figure [5]d has a characteristic timescale 
that is apparently comparable to the exposure time, and 
thus the structure function amplitude is flat over the 1-3 
hour time separations that are sampled. In this case, the 
structure function amplitude is consistent between the 
SDSS and Hydra samples. 

The increase in variability amplitude with time seen 
in our low -variability (high activity) sample is reminis- 
cent of the lLee et al.l (|2010[ ) result that there was a much 
larger number of Ha emission events that lasted at least 
30 minutes compared to the number of shorter events 
they observed. Their sample was primarily composed of 
very active stars which is comparable to this subset of 
our data. 

Following the discussion in §3.2, we interpret the 
timescale results as follows. The low-variability (high 
activity) stars are covered with spots, so that small vari- 
ations are not easily visible against the high background. 
In order to see a noticeable variation, a rather significant 
event must occur. Since higher-energy events are known 
to las t longer (e.g. iLacv et al]|1976l : iHawlev &: PettersenI 
|1991|) , this leads to a longer timescale for events that are 
strong enough to be detected above the persistent activ- 
ity. For the high-variability (low activity) stars, small 
events will cause a noticeable change, and these occur 
on short timescales, below the ~15 minute limit of our 
sampling cadence. 

4. SUMMARY 

We used two spectroscopic samples of active M dwarfs 
to examine Ha emission variability on timescales from 



minutes to weeks. 

We found that the Ha variability measured us- 
ing the Rew metric remains relatively constant as a 
function of spectral type, in contrast to previous re- 
sults. However, using our fractional variability met- 
ric (cr£;vi/ffyi/(EWHA)), we showed that the higher- 
variability stars have relatively low activity strength 
(LHa/L&o;)j conversely the low- variability stars have 
high activity strength. This naturally leads to an appar- 
ent dependence on spectral type because the later type 
stars typically have lower activity strength. We specu- 
lated that the physical reason for the higher variability 
in the low activity stars is that small changes in Ha emis- 
sion, e.g. due to small-scale flaring or the appearance of 
a new active region on the stellar surface, would have a 
relatively larger effect on the less active stars. 

We investigated the timescales of Ha emission variabil- 
ity using structure functions. The observed Ha varia- 
tions for low- variability stars (cr_EH'ffA/(EWHA) < 0.16) 
occur on a timescale longer than an hour, while the high- 
variability M dwarfs exhibit a timescale that is shorter 
than our sampling time of ^15 minutes. Neither the low 
nor high variability samples show significant changes in 
the structure function on timescales longer than a day. 
We speculated that the low- variability sample, which is 
mainly composed of the higher activity strength stars, 
requires more energetic and hence longer-lived emission 
events in order to be detected above the strong persis- 
tent surface activity, leading to a longer characteristic 
variability timescale. 

Better time resolution in spectroscopic M dwarf mon- 
itoring, as well as a wider and better-sampled range of 
time separations, will allow for a more accurate deter- 
mination of the timescales of Ha variability in low-mass 
dwarfs and may lead to a better understanding of the 
physical mechanisms that cause these changes. 
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